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M. Kugler,i’2,[^G. Brandi,J. Waizner,^ M. Janoschek,^ R. Georgii,k2 
A. Bauer,^ K. Seemann,^’^ A. Rosch,^ C. Pfleiderer,^ P. Boni,^ and M. Garst^ 

^Physik-Department E21, Teehnisehe Universitdt Miinehen, 85748 Garehing, Germany 
^Heinz Maier-Leibnitz Zentrum (MLZ), Teehnisehe Universitdt Miinehen, 8574-8 Garehing, Germany 
^Institut fiir Theoretisehe Physik, Universitdt zu Kdln, Ziilpieher Str. 77a, 50937 Kdln, Germany 
Gondensed Matter and Magnet Seienee, Los Alamos National Laboratory, Los Alamos, New Mexieo 87545, USA 

(Dated: August 31, 2015) 

A magnetic helix realizes a one-dimensional magnetic crystal with a period given by the pitch 
length Xh. Its spin-wave excitations - the helimagnons - experience Bragg scattering off this pe¬ 
riodicity leading to gaps in the spectrum that inhibit their propagation along the pitch direction. 

Using high-resolution inelastic neutron scattering the resulting band structure of helimagnons was 
resolved by preparing a single crystal of MnSi in a single magnetic-helix domain. At least five he- 
limagnon bands could be identified that cover the crossover from flat bands at low energies with 
helimagnons basically localized along the pitch direction to dispersing bands at higher energies. In 
the low-energy limit, we find the helimagnon spectrum to be determined by a universal, parameter- 
free theory. Taking into account corrections to this low-energy theory, quantitative agreement is 
obtained in the entire energy range studied with the help of a single fitting parameter. 

PACS numbers: 75.30.Ds, 78.70.Nx, 71.27.-ha, 71.70.Gm 


The weak spin-orbit Dzyaloshinskii-Moriya interac¬ 
tion, in the cubic chiral magnets energetically favours 
spatial modulations of the magnetization. This gives rise 
to magnetic crystalline phases with unit cells that are 
incommensurate with and much larger than the atomic 
lattice spacing. Most prevalent is the magnetic helix, 
a one-dimensional magnetic crystal, with a large pitch 
Xh = ^nlkh proportional to the ratio J/D where J is the 
magnetic exchange [T]. For a small range of finite mag¬ 
netic fields, a two-dimensional magnetic crystal is also 
stabilized close to the critical temperature m- It can 
be identified as a lattice of magnetic skyrmions whose 
non-trivial topology is at the origin of various interesting 
phenomena [3] like, for example, a topological Hall effect 
gug and an emergent electrodynamics EllTl. Interest¬ 
ingly, the phase transition from the paramagnetic to the 
magnetically ordered phases at small fields corresponds 
to a weak crystallization process [8] and is driven first- 
order by strongly correlated chiral paramagnons [9HII]. 

The spin-wave excitations of these magnetic crystals 
possess a band structure uJn,ci with band index n that, 
according to Bloch’s theorem, reflects the periodicity of 
the magnetic order. For the magnetic helix with a pitch 
vector k/j,, the dispersion is periodic, = ^n,q+mk^ 

with m G Z, along the direction in momentum space 
singled out by k/,,. In contrast to commensurate anti- 
ferromagnets, however, the size ky = |k/i| of the result¬ 
ing magnetic Brillouin zone is small; for MnSi at lowest 
temperatures Xy = 180A and ky = 0.035A“^. Impor¬ 
tantly, this ensures on the one hand that the dispersion, 
ujn,q^ of the magnons is universal in the sense that it is 
captured by an effective continuum theory and is deter¬ 
mined by only a few parameters. On the other hand, 
a high resolution of momenta is required in order to re¬ 


solve the band structure experimentally with the help 
of inelastic neutron scattering. In a first experiment on 
MnSi, Janoschek et al m succeeded to acquire scatter¬ 
ing spectra for the helimagnons and also described them 
theoretically but neglecting dipolar interactions Haul]. 
In MnSi at zero magnetic field, the pitch vector of the 
helix, k/j,, aligns with one of the eight equivalent crys¬ 
tallographic (111) directions giving rise to four magnetic 
domains each of which is determined up to a phase. At 
that time, the experiment was performed on a single crys¬ 
tal containing multiple domains so that magnon excita¬ 
tion branches of all four domains were simultaneously 
excited. This resulted in broad total spectra rendering 
the identification of individual magnon modes impossible, 
see Fig. Bf)- Individual modes, however, have already 
been identified with the help of magnetic resonance ex¬ 
periments on MnSi, Cu 20 Se 03 , and Feo. 8 Coo. 2 Si dSHIZl 
that probe certain branches but only at the center of the 
Brillouin zone, i.e., at frequencies co’n,q=o- Schwarze et al. 
m were able to explain these resonances quantitatively 
after taking dipolar interactions into account. 

In the present study, the band structure of the heli¬ 
magnons could be resolved by first preparing a crystal 
of MnSi in a single magnetic-helix domain by applying a 
small field. A finite magnetic field H competes with the 
crystalline anisotropies to align the pitch vector, k/^, of 
the helix. For sufficiently large H > Hd , the pitch vector 
is parallel to H as this allows the magnetization to cant 
towards the field gaining magnetic Zeeman energy giving 
rise to a single magnetic-helix domain. The magnetiza¬ 
tion M = mn with = I is then governed by the free 
energy density T = T^mag + ‘^dipolar where the first and 
second term account for the short-range magnetic and 
the dipolar interactions, respectively. In the low-energy 
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FIG. 1: (a)-(e) Universal helimagnon spectrum for momenta q\\ along the pitch direction, k/^, in the repeated zone seheme (thin 
lines) at zero field and Xcon = 0-34 for various values of perpendicular momentum q±. The thickness of the colored lines indicate 
the theoretically expected intensity for inelastic neutron scattering on a domain with k/^ || [111] with respect to a nuclear (110) 
reflection and for || [112]. The colors red and green distinguish between dispersive and practically flat helimagnon bands, 
respectively. The energy units on the right axes correspond to MnSi at 20 K for which = 0.062 meV. (f) Weights 

of inelastic contributions arising from an equal population of all domains (color coding as in (g)) for q±_ — 3kh as a function of 
qy = ^||k/i with respect to k/^ || [111] (red domain), (g) Scattering plane around the nuclear (110) peak with four domains that 
are degenerate at zero field; the open circles correspond to out-of plane pitch wavevectors projected onto the plane. 


limit the former reduces to J-mag ~ J~o with the standard 
model for chiral magnets [1] 



(Vifij)^ + 2khon{V x n) 


— iiomhia. 


( 1 ) 


Minimization of Tn + ‘^dipolar yields the conical helix 
Uq = (sin^cos(/c/,, 2 :), sin^sin(/c/i^), cos^) for an applied 
field H along the 2 ;-axis with kh = kho- The angle 0 
parametrizes the homogeneous part of the magnetization 
that depends linearly on the magnetic field mcosO = 
xi^onHint where the internal field = H/{1 + N^Xmn) 
with demagnetization factor Nz arising from /bipolar and 
the susceptibility being related to 

the stiffness density The transition to the field- 

polarized state then occurs at the second critical field 
Hlf = m/x“n = Pskl/ipom). 

Plugging the standard parametrization fi = 
63^1 — 2|'0P + '06+ + 0*6“ with 6+ = (ei ± 162 ) 1^/2 
and where 63 = no into the Landau-Lifshitz 

equation and expanding in lowest order in 0 one finds 
that for intermediate fields Hd < H < i^c 2 , the magnon 
wave-function 0 is governed by a bosonic Bogoliubov- 
deGennes equation ihr^dt^ = 1-L^ for the spinor 

= (0,0*). The Hamiltonian 71 = Tio 'Hdipoiar 
consists of two parts that derive from Tq and Jyiipoiar, 
respectively. The former contribution reads [12 HU 


no = v 


-IV^ - t2T^khn±{r)\/ + -(1 


( 2 ) 


where and are Pauli matrices, and the stiffness is 
V = giJ.Bps/'m = p ^ 2 for MnSi. 

The helimagnons are subject to an effective vector po¬ 
tential n^(r) = {sinO cos{khz),sm0 sm{khz),0) that is 
periodic and is responsible for the formation of bands. 

In the low-energy limit the Hamiltonian depends only 
on three parameters: (i) the critical field 


pitch vector k/j,, and (in) the numerical value Xcon = 0-34 
for MnSi that effectively measures the strength of the 
dipolar interactions di. All these parameters are known 
from independent measurements resulting in a universal 
prediction for the helimagnon spectrum that we discuss 
in the following. 

For momenta q = q\\i^h strictly longitudinal to the 
pitch vector, the helimagnons decouple from the periodic 
potential in Eq. The helimagnon dispersion can then 
be obtained in closed form and reads in the extended zone 
seheme^ see Fig. [^a), 

fiw(q = qi\kh) = V\q\i\yjqf^ + (1 + ^ (3) 

where sin^ ^ = 1 — (i7int/^c2^)^ and the dependence on 
Xcon attributed to the dipolar interaction. For large 
momenta \c[± \ kh perpendicular to the pitch, q^k^ = 
0, on the other hand, the dipolar interaction as well as 
the last term in Eq. (§ can be neglected and the wave 
equation reduces to a Schrodinger equation 




(4) 


V 


qj_- + 2khq± sin 9 cos{khZ - a) 'ipoj,ci_L (z) 


where {qx^Qy) = q±{cosa,sma). It effectively describes 
a particle in the presence of a periodic cosine potential 
along the 2 :-axis, i.e., it is the Mathieu equation. Inter¬ 
estingly, the strength of the potential can be tuned by q±_ 
with the concomitant reconstruction of the spectrum as 
illustrated in Eigs. Bb) -(e). A finite q^ activates Bragg 
reflections that open gaps at the Bragg planes. As the 
potential only contains the two primary Eourier compo¬ 
nents, the size of these band gaps however strongly de¬ 
creases as a function of the band index n = 0,1,2,... 
|18| . Eor large q± ^ {n l)‘^kh/ sin^ the potential for a 
given band n eventually becomes strong and basically lo¬ 
calizes the helimagnon along the pitch direction. In this 
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tight-binding limit, the bands are practically flat, i.e., 
independent of gy, and the spectrum is given by 


^^n,q ^ ^ 


q]_ - 2khqi. sin 0 (5) 

_ / \ 

+ 2kh\/kuqi. sin0 \^n+ -j - + n) 


up to corrections of order 0{Vk\). The second and third 
term derive from an expansion of the cosine potential up 
to quadratic order giving rise to a harmonic oscillator 
spectrum; the last term is attributed to the anharmonic- 
ity of this potential. 

In order to verify the helimagnon theory, we have 
conducted inelastic neutron scattering experiments at 
the cold neutron triple-axis spectrometer MIRA-2 m 
at the neutron source FRM-II in Garching using neu¬ 
trons with fixed incident energies Ei = 3.19, 4.06 and 
5.04 meV. MIRA-2 is particularly well suited for the ex¬ 
periments due to the intrinsic excellent momentum reso¬ 
lution, which is important for resolving spin waves close 
to magnetic satellite peaks. The collimations before and 
after the sample were 30’ resulting in energy resolutions 
AE = 76, 118 and 161 /ieV, respectively. Higher order 
neutrons were removed by a cooled Be-filter. Two large 
single crystals of MnSi were used in the present study 
with a volume of approximately 8 cm^ and a small mo- 
saicity of 10’. 

For the measurement of the helimagnon spectrum 
three experimental setups were realized, see Fig. [la)-(c). 
A magnetic field H > Hd was applied in all setups to 
align the pitch vector and, more importantly, to prepare 
a single magnetic-helix domain. In setup 1, Fig. |^a), 
the magnetic field /io^int ~ 0-07 T was within the scat¬ 
tering plane and pointing along the [111] direction mea¬ 
sured with respect to the nuclear (110) Bragg peak. This 
configuration is ideal to obtain an excellent energy reso¬ 
lution for the dispersive branches along gy with g^ = 0. 
In setup 2, Fig. |^b), a magnetic field /io^int ~ 0-13 T 
was applied perpendicular to the scattering plane along 
[112] aligning the pitch with [112], and the helimagnons 
were measured for momenta along [111]. In setup 
3, Fig. ^c), measurements were performed on the sec¬ 
ond crystal as a consistency check with a magnetic field 
/^o^int ~ 0-13 T applied perpendicular to the scattering 
plane along the [110] direction orienting the pitch along 
[110], which allowed us to probe the bands for along 
[001]. The data acquired along [001] and [111] turned out 
to be identical. 

Most measurements were performed at 20 K where the 
intensity of the magnetic Bragg peaks is high and the 
Bose factor is large yielding a high inelastic intensity. 
As an example. Fig. [^d) shows a sequence of scans ob¬ 
tained with setup 1 with momenta longitudinal to k/^, 
i.e., q± = 0. According to Fig. [^a) weights from three 
branches are expected and can indeed be identified for 
gy = 2kh- The high-energy mode near /kj ~ 0.5 meV is 


(a) Setup 1 


(b) Setup 2 


(c) Setup 3 
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(f) Setup 3 



FIG. 2: Panels (a)-(c) show the scattering planes for the 
setups 1, 2, and 3, respectively. The magnetic field (blue 
color code) defines the orientation of k/^ resulting in magnetic 
satellite peaks (green squares); for panel (b) and (c) they are 
located above and below the scattering plane. Constant-Q 
scans have been performed at the points in reciprocal space 
marked by red circles. Panels (d)-(f) show examples of scans 
obtained with the three setups at 20 K, where the full lines 
are fits of multi-Gaussian profiles to the data, (d) Three heli¬ 
magnon branches are measured with setup 1 (incident energies 
Ei = 3.19 and 4.06 meV), while (e) setup 2 (Ei = 4.06 meV) 
and (f) setup 3 {Ei = 5.04 meV) clearly resolve the first two 
and five bands, respectively. The elastic peaks appear due to 
the Q-independent incoherent scattering. Note that the scans 
in (d) are shifted by 40 counts with respect to each other for 
clarity. 


best seen for incident energy Ei = 4.06 meV {O symbols) 
and the other two are resolved for Ei = 3.19 meV (O 
symbols). While the low-energy mode at fuj ^ 0.1 meV 
is partly hidden by the incoherent scattering, the mode 
with intermediate energy at hw ~ 0.25 meV is most 
prominent. Both shift to larger energies with increas¬ 
ing gy. Fig. j^e) shows the g^-dependence of the first 
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(a) Setup 2 & 3, T = 20 K 


(b) Setup 1, T = 20 K 


(c) Temperature Dependence 



FIG. 3: Comparison between theory and experiment (symbols), (a) & (b) Helimagnon spectrum as a function of q±_ and q\\ 
with nominally q\\ — 0 and = 0, respectively. The universal spectrum (dotted lines) agrees with the data at low energies; the 
correction due to Eq. 0 is necessary (solid lines) for quantitative agreement at higher energies. The instrumental resolution 
slightly shifts the lower and center branch upwards (black dashed line in panel (b)), see text, (c) T-dependence of the two 
helimagnon modes lowest in energy measured at certain fixed momenta. 


two bands which are accessible with setup 2 at gy = 0. 
With increasing q±_ the energy of the bands increases as 
expected. Finally using setup 3, we succeeded to probe 
the annihilation of helimagnons that belong to at least 
five different bands for qx = —3kh{0,0, 1), see Fig. |^f). 

For a quantitative analysis, the data was fitted to a 
multi-Gaussian spectrum after subtracting a background 
of 1.0 cts/min and 1.5 cts/min for setup 1 and se¬ 
tups 2 & 3, respectively. The latter was determined by 
repeating most measurements at 3 K, where the inelas¬ 
tic magnetic scattering is weak. The fits are shown as 
solid lines in Figs.[^d)-(f). Note that the observed width 
varies for different peaks because the resolution of the 
spectrometer depends on the energy transfer Huj. The 
peak positions obtained from these fits are summarized 
in Fig. 

Fig. I^a) and (b) compares the experimental data 
at 20 K to the prediction of the universal helimagnon 
spectrum (dotted lines) using ~ 0.036A“^ [20] and 
^ 0.53 T [2T], that determines the energy scale 
= 0.062 meV at 20 K. While there is perfect 
agreement at low energies, deviations become substantial 
at higher energies with the data assuming systematically 
lower values than predicted. It turns out that these devi¬ 
ations can be attributed to corrections to the low-energy 
theory 0 caused by higher-order contributions in the 
gradient expansion. Representatively, we consider the 
following correction [22| 

<5Jmag = f (6) 

where ^ is a dimensionless number that is expected to be 
small and of order A ^ 0{{akhY) with ^ 24 being 
the volume of a formula unit in MnSi so that {ak^Y ~ 
0.01. For momenta g ~ 3/c/,,, for example, the corrections 
arising from Eq. 0 are expected to be 3‘^{akh)‘^ ^ 10% 


as compared to the first term in Eq. Q, which is in fact 
comparable to the deviations observed in Fig. For a 
quantitative comparison in the entire investigated energy 
range, we took the modification of the helimagnon spec¬ 
trum due to Eq. 0 into account providing us with a 
single fit parameter A. The solid lines in Fig. [^a)- (c) 
show a best-fit yielding Am = —0.0073 ± 0.0004. 

The bands in Fig.j^a) are mostly in the tight-binding 
limit and thus practically independent of gy. As a result, 
the measured peak positions are rather insensitive to the 
vertical momentum resolution of MIRA (Ag « ±0.6kh)- 
The flatness of the lowest bands has been also explic¬ 
itly checked by comparing the two scans at finite q± in 
Fig. [^a) using setup 1 (not shown). In contrast, the dis¬ 
persive bands in Fig. [^b) have a substantial q± depen¬ 
dence which is sampled by the instrumental momentum 
resolution. The resulting admixture of spectral weight 
with finite q± into the nominal q± = 0 spectrum leads to 
a slight upward shift of the expected peak positions as 
indicated by the black dashed line in panel Fig.|^b) |23| . 

In order to investigate the temperature dependence, 
the helimagnons were measured at gy = 2 .1kh with setup 
1 and at q± = 2.5kh with setup 2 at four and six dif¬ 
ferent temperatures, respectively. Fig. |^c) shows the T- 
dependence of the peak positions for the two helimagnon 
modes with lowest energies. The T-dependence mainly 
derives from the magnitude of the local magnetization 
m(T) that is reflected in the critical field enter¬ 

ing the stiffness V in Eq. 0 . The theoretically expected 
peak positions resulting from the values of ob¬ 

tained from ac susceptibility measurements m as well 
as from the weakly T-dependent pitch length kh{T) [20| 
are shown as solid lines in Fig. |^c) in good agreement 
with experiment. 

In the present work, we experimentally verified the the¬ 
ory for helimagnons and established the emergence of flat 
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helimagnon bands at finite momenta. The latter might 
be technologically exploited, for example, in the design of 
magnon waveguides in magnonic applications [24]. Our 
findings set the stage for further investigations of heli- 
magnons at finite magnetic field and pressure and, in 
particular, their influence on the stabilization of various 
phases EEHII]- It has been suggested that the renormal¬ 
izations due to magnons are crucial for the stabilization 
of the skyrmion crystal [5|. They might be also impor¬ 
tant for the formation of the enigmatic non-Fermi liquid 
phase [23123 observed in MnSi at high pressures [27l|28|, 
a long-standing puzzle of this material that still remains 
unsolved. 
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